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Abstract 

Using data collected in the DELPHI detector at LEP-1, measurements of the in- 
clusive r branching ratios for decay modes containing one, three, or five charged 
particles have been performed, giving the following results: 

B x = B(t~ -»■ (particle)" > Ovr > 0K°v T (v)) = (85.316 ± 0.093 ± 0.049)%; 
B 3 = B{t~ -> 2h-h + >0ir°>0K°v T ) = (14.569 ± 0.093 ± 0.048)%; 

B 5 = B{t- ^3h-2h + >0n >0K°u T ) = (0.115 ± 0.013 ± 0.006)%, 

where h is either a charged 7r or K meson. The first quoted uncertainties are 
statistical and the second systematic. 
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1 Introduction 

Historically, there has been an inconsistency in the measured r exclusive branching 
ratios: the measurements of exclusive branching ratios to decay modes containing one 
charged particle did not sum up to the inclusive branching ratio for a charged multi- 
plicity of one. A number of analyses have attempted and succeeded in resolving this 
question HHHI^a] . However the current world average values from direct measurements of 
the topological branching ratios in the Particle Data Group listings || have uncertainties 
which are significantly larger than the values obtained through combined fits to all the r 
decay data. In the case of the 1-prong branching ratio (T2 in ||), the difference between 
the average and the fit value is more than two standard deviations. 

This paper presents a dedicated simultaneous measurement of the decay rates of the 
r lepton to different final states as a function of the charged particle multiplicity. 

The relevant components of the DELPHI detector and the data-set are described in 
Section 0. The definition of the measured branching ratios and the method used to 
derive them are introduced in Section |3|. In Section |] the preselection of the sample of 
e + e~ — ► t + t~ events is described. Section [5] contains a description of the reconstruction of 
charged particles, as well as the reconstruction algorithms for photons and hadrons which 
interact with the detector material before the tracking subdetectors. The r decays are 
classified according to their charged particle multiplicity. This is described in Section | 
The results and systematic studies are presented in Sections BTTI and FT2I respectively. 



2 Experimental apparatus and data sample 

The DELPHI detector and its performance are described in detail elsewhere H0. 
The subdetector units particularly relevant for this analysis are summarised here. These 
detector components covered the full solid angle considered in the analysis except where 
specified, and sat in a 1.2 Tesla solenoidal magnetic field parallel to the z-axis[]. 

The reconstruction of charged particles in the barrel region of DELPHI used a com- 
bination of the measurements in four different cylindrical subdetectors: a silicon Vertex 
Detector (VD); the Inner Detector (ID), consisting of a jet chamber tracker and a wire 
chamber used for trigger purposes; the Time Projection Chamber (TPC); the Outer 
Detector (OD). 

The VD had three layers of silicon micro-strip modules, at radii of 6.3, 9.0 and 11.0 cm 
from the beam axis. The space point precision was about 8 fim in r-0. For data from 
1993 onwards, a measurement of r-z was produced in the outermost and innermost layers. 
This had a precision of about 15 /im. The-two track resolution was 100 /im in r-0 and 
200 /im in r-z. 

The ID had an inner radius of 12 cm and an outer radius of 28 cm. The inner jet 
chamber part lay at radii below 22 cm. It had a two-track resolution in r-0 of 1 mm and 
a precision in r-<fi of 50 /im. The outer triggering part was a five-layer wire chamber. This 
was replaced for the 1995 data with a straw detector containing less material. 

The TPC, extending from 30 cm to 122 cm in radius, was the main detector for charged 
particle reconstruction. The main track reconstruction information was provided by 16 
concentric circles of pads which supplied up to 16 three-dimensional space points on a 
track. In addition, ionisation information was extracted from up to 192 wires. This was 

1 In the DELPHI reference frame the origin is at the centre of the detector, coincident with the interaction region, 
the z-axis is parallel to the e _ beam, the x-axis points horizontally towards the centre of the LEP ring and the y-axis is 
vertically upwards. The co-ordinates r,</>,z form a cylindrical coordinate system, while 8 is the polar angle with respect to 
the z-axis. 
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used for particle identification purposes. Every 60° in (j> there was a boundary region 
between read-out sectors about 1° wide which had no instrumentation. At cos# = there 
was a cathode plane which caused a reduced tracking efficiency in the polar angle range 
| cos#| < 0.035. The TPC had a two-track resolution of about 3 mm in r-0 and 1.5 cm 
in z. 

The OD, with 5 layers of drift cells at a radius of 2 m from the beam axis, was 
important for the momentum determination of energetic particles. 

The principal device for electron and photon identification was the High density Pro- 
jection Chamber (HPC), located outside the OD, which allowed full reconstruction of 
the longitudinal and transverse shower components. It covered the polar angle region 
| cos#| < 0.75. In the forward region, 0.800 < | cos#| < 0.985, the calorimetry was per- 
formed by a lead-glass array. For very small polar angles, the calorimetry was performed 
by the luminosity monitors, the SAT in 1992 and 1993, and the STIC in 1994 and 1995. 
In the region 0.75 < | cos#| <0.8, there was no electromagnetic calorimeter. 

Outside the magnet solenoid lay the hadron calorimeter and muon chambers, which 
were used for hadron and muon identification. 

The data were collected by the DELPHI detector from the LEP electron-positron 
collider, in the years 1992 through 1995, at centre-of-mass energies y/s of the e + e~ system 
between 89 and 93 GeV, on or near to the Z resonance. It was required that the VD, 
ID, TPC and HPC were in good operational condition for the data sample analysed. 
The integrated luminosity of the data sample was 135 pb _1 of which about 100 pb _1 was 
taken at yfs of 91.3 GeV, close to the maximum of the Z boson resonance production 
cross-section. 

The selection procedures were studied with samples of simulated events which had 
been passed through a detailed simulation of the detector response J7J and reconstructed 
with the same program as the real data. Separate samples were produced correspond- 
ing to the detector configurations in different years. The Monte Carlo event generators 
used included: KORALZ 4.0 f§ for e+e" -> t+t~ events; DYMU3 § for e+e" -> 
events; BABAMC |jj and BHWIDE Q for e+e" -> e+e" events; PYTHIA 5.7 H 



for e + e~ — > qq events. Four-fermion final states were produced using two different gen- 
erators. BDK |13| was used for reactions with four leptons in the final state. This 
included two-photon events where one or two e + or e~ were not observed in the detector. 
TWOGAM [14] was used to generate e + e~ — > (e + e~)qq events. The KORALZ 4.0 gen- 
erator incorporated the TAUOLA 2.5 |H| package for modelling the r decays. A value 
of 1.2% was used for the branching ratio of the Dalitz decay 7r° — > je + e~ . 



3 Method 

From a high purity sample of Z — ► t + t~ events, r decays were identified with charged 
multiplicity one, three or five, corresponding to the topological branching ratios de- 
fined as: 

B x = B(r~ -> (particle)" >07r >0ir°z/ r (z/)); 
B 3 = B(t~ -> 2h-h + >07i°>0K°v T ); 
B 5 = B(t~ -> 3h-2h + >0n°>0K°v T ), 

where h is either a charged ir or K meson. In this definition, the Kg meson was treated 
as a neutral particle, even when it decayed to 7r + 7r _ before the tracking detectors. In 
r decays only about 5% of K$ mesons decay before the VD. The h~ K° (neutral s)v T 
component was therefore treated as signal for the 1-prong r decays and background for 
the 3-prong r decays. Although a tt° meson which decayed to the 7e + e~ final state 
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produced two charged particles observed in the detector, it was also treated as a neutral 
particle, following the r branching ratio definition given by the Particle Data Group ||. 
Thus the branching ratios B\, B 3 and B 5 correspond to the parameters T 2 , r 55 and Tioi 
respectively in 

No attempt was made to measure r decays to more than five charged particles and it 
was assumed that such decays had a negligible branching ratio. A 90% CL upper limit 
on the inclusive branching ratio for the r decay to seven charged particles has been set 
at 2.4 x 10~ 6 by the CLEO experiment |fj. 

At LEP it is possible to separate cleanly t + t~ events from other final states, permitting 
the efficient selection of a high purity r sample. One can measure the branching ratio 
B(t — > X) for the decay of the r to a final state X using the expression 

where Nx is the number of identified decays of type X found in the sample of N T r-decay 
candidates, preselected with efficiency e T with a background fraction of b T . ex is the 
total (preselection x identification) efficiency for selecting the decay mode r — > X, with 
a background fraction of bx, including background from other r decay modes. In a t + t~ 
event selection without specific requirements on either of the two r candidates in the 
event, e r will be identical to the t + t~ event selection efficiency. 

In the case where several branching ratios are measured simultaneously, candidate r 
decays can be classified as a function of the detector components used in the charged 
particle track reconstruction. These classes (detailed in Section |5.4|) typically have dif- 



ferent signal to background ratios for different types of decays and thus contain different 
levels of information. To obtain the branching ratios, a fit can then be performed to the 
predicted number N i)Pred of decays in a class i: 

Ni tPre d = N T {ei t \B\ + e it3 B 3 + e it5 B 5 ) + N i:bkg , (2) 

where 6i jm is the probability of a r decay of true multiplicity m being attributed to class 
i and N^g is the estimated background in class i due to non-r + r~ events. This can 
be extended to a simultaneous fit for both r decays in a candidate t + t~ event. The 
predicted number N it j tPre d of t + t~ events with one r decay in class i and the other in 
class j is 

^i,j,pred X TT ^ ] 6^ B m B n -\- Nijfikg 5 (3) 
m,n 

where Nij^kg is the estimated number of non-r + r~ background events in class i, j and e^ n 
is the probability of an event with decays of true multiplicity m and n being attributed 
to class N TT is the number of r + r~ events and can be left free and obtained from a 
fit to the different classes. This method takes into account correlations between the two 
observed r decays in an event and has the advantage that non-r + r~ backgrounds tend 
to populate event classes which have correlations between the two candidate r decays 
(e.g. dimuons in the 1- versus- 1 topology but not in l-versus-3) and thus exhibits less 
sensitivity to backgrounds. It was used as the principal method in the following analysis. 



4 e + e — > t t preselection 

In e + e~ — » Z — ► r + r~ events at y/s = M z , the r + and r~ are produced back-to-back, 
ignoring radiative effects. Each r decays, producing one, three, or more charged particles 
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in addition to one or two neutrinos and, possibly, neutral mesons. All particles apart from 
the neutrinos can be detected in DELPHI. A t + t~ event is thus characterised by two low 
multiplicity jets which appear approximately back-to-back in the laboratory. The loss of 
the neutrinos implies that not all the energy in the event is seen in the detector. 

Except where explicitly stated, all quantities calculated using charged particles used 
only charged particle tracks which had hits in at least the TPC or OD, to ensure good 
momentum reconstruction from the long track length, and had an impact parameter with 
respect to the centre of the interaction region of less than 1.5 cm in r-0 and 4.5 cm in z. 
The basic e + e~ — > t + t~ preselection described below has some differences compared 
with that used previously in studies of the r leptonic branching ratios jl7| . There are 
slightly modified cuts to reduce non-r backgrounds and the addition of subsamples of 
events containing r decays with a nuclear reinteraction in the detector material and r 
decays where one charged particle was reconstructed only in the VD and ID tracking 
subdetectors. 



4.1 Preselection criteria 

To ensure that the r decay products were in the region of DELPHI corresponding to 
the HPC acceptance, the thrust axis of the event, calculated using only charged particle 
tracks, was required to lie in the polar angle region defined by | cos#| < 0.732. The event 
was split into two hemispheres, each associated to a candidate r decay, by the plane 
perpendicular to the thrust axis and containing the point at the centre of the interaction 
region. Each hemisphere had to contain at least one charged particle. It was required 
that at least one charged particle lie in the polar angle region defined by | cos#| > 0.035 
in order to reduce the effects of the TPC tracking inefficiency near cosO = 0. 

Most hadronic Z decays were rejected by requiring that the event contain at most 
eight charged particle tracks with hits in either the TPC or OD and satisfying the impact 
parameter cuts described above. 

Four-fermion events were suppressed by requiring that the event isolation angle be 
greater than 160°. This was defined as the minimum angle between any pair of charged 
particles in opposite r decay hemispheres. 

Backgrounds from n + fi~ and e + e~ final states and cosmic rays were reduced by re- 
quiring that the isolation angle be less than 179.5° for events with only two charged 
particles. 

The and e + e~ contamination was reduced further by requiring that both p Ta d = 

f^i- + 1/2 and E rad = {^k + -ffe) 1/2 be less than unity. The variables p\ and p 2 are the 
momenta of the highest momentum charged particles in hemispheres 1 and 2 respectively. 
The quantity p[ was obtained from the formula 

pi = y/ssme 2 /(sme 1 + sinfl 2 + | sin(0i + 2 )|) , (4) 

and p' 2 by analogy with the indices 1 and 2 interchanged. The angles B% and 9 2 are the po- 
lar angles of the highest momentum charged particle in hemispheres 1 and 2 respectively. 
The variables E% and E 2 are the total electromagnetic energies deposited in cones of half- 
angle 30° about the momentum vectors p\ and p 2 respectively, while Ej = cp'p for j = 1,2. 
These definitions of p ra( i and E Ta( i are identical to those used for the e + e~ — > t + t~ mea- 
surements at LEP-2 JTB|. They take into account more correctly the radiative events 
where the radiative photon is close to the beam axis than the definition used previously 
in LEP-1 analyses [|i~7| , |T9l . In the limit of no radiation the new and old definitions are 
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Figure 1: Distributions of a) p ra d and b) E rad variables used in t + t~ selection. The 
dots are data, the line is simulation. The shaded area is non-r + r~ background. In both 
cases the selected sample lies in the region below a cut at 1.0. The spike in the first bin 
of plot (b ) is due to events where both r 's decay to either a muon or a charged hadron 
which does not interact in the UPC. This bin has been scaled down by a factor two for 
presentational purposes. 

equivalent. The distributions of p ra d and E ra d are shown in Fig. pi and Fig. |l|b respec- 
tively, after all cuts have been imposed except the one on the displayed variable. There 
are some small data/simulation discrepancies in the momentum resolution, visible in the 
dimuon peak at p ra d ~ 1-4. This had a negligible effect on the analysis. 

Much of the remaining background from the dileptonic channels came from events 
containing hard radiation lying far from the beam. These events should lie in a plane. 
Where two charged particles and a photon were visible in the detector, such events were 
removed by requiring that the sum of the angles between the three particles was greater 
than 359.8°. 

A further reduction in four-fermion background was achieved by requiring a minimum 
visible energy in the event. The visible energy E vis is the sum of the energies of all 
charged particles and the electromagnetic energies of all neutral particles, neglecting 
energy deposits recorded by the very forward calorimeters (SAT and STIC) at angles of 
less than 12° from the beam axis. Events were accepted if E V i S was greater than 0.09 x ^/s. 
For events with only two charged particles, an additional condition was applied that the 
vectorial sum of the charged particle momenta in the r-0 plane be greater than 0.4 GeV/c. 
Four-fermion events typically have very low values of of this quantity compared with t + t~ 
events. 

Most cosmic rays were rejected by the upper cut on isolation angle. Further rejection 
was carried out by requiring at least one charged particle in the event have an impact 
parameter with respect to the interaction region of less than 0.3 cm in the r-0 plane and 
that in the r decay hemisphere opposite to this particle at least one charged particle have 
an impact parameter in the r-0 plane of less than 1.5 cm. It was also required that both 
event hemispheres have a charged particle track whose perigee point lay within 4.5 cm of 
the interaction region in z. 
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Figure 2: The output neuron distribution for the neural network used for the final qq 
rejection. The dots are data, the line is simulation and the shaded region shows the 
simulated qq contribution. The selected events lie in the region with output neuron values 
above 0.05. 



Additional r + r~ events were selected where one of the r's produced a 1-prong decay 
and the other r decay produced one charged hadron which interacted with the detector 
material inside the TPC and where the interaction was reconstructed by the nuclear 
reinteraction reconstruction algorithm described in Section |5.3| . The cuts in p rac i and 
E ra d were not applied. This increased the selection efficiency by about 1% with a relative 
background of 0.8%. 

A further class of events was selected where one of the r decays left one charged 
particle track with hits in only the VD and the ID due to inefficiencies in the TPC and 
OD. Many of these tracks contained no polar angle information and had poor momentum 
resolution, so they were excluded from the calculation of event quantities such as the 
thrust or p ra d- The p ra d requirement was replaced with a cut on the momentum of 
the highest momentum charged particle in the opposite hemisphere, requiring it to be 
less than 90% of the beam momentum. This cut was tightened to 70% if there was an 
identified muon in the hemisphere. This subsample of events contributed an extra 4% of 
events with a relative background of 5%. 

The total efficiency, including all classes of events, was estimated from simulation to 
be (51.74 ± 0.04)%. This corresponds to an efficiency of approximately 85% within the 
angular acceptance cuts. 

4.2 Backgrounds 

Knowledge of the background from qq events is of particular importance in this analysis 
as it constitutes a serious background in the higher multiplicity topologies. Uncertainties 
in the low multiplicity backgrounds have a reduced effect on the measured branching ratios 
due to the large B\ value and the resultant cancellation in Eqn. |l|. After the cuts described 
above, the qq background estimated from simulation is (0.80 ± 0.05)%, where the error 
is purely statistical. To reduce systematic effects due to qq background, a further level 
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Figure 3: Distributions of isolation angle (in degrees) in candidate t + t~ events: a) events 
containing two identified muons; b) events containing at least one tagged electron; c) all 
events in the t + t~ sample. In all cases data are the dots and simulation is the line. In 
(a) and (b) the shaded area is the four-fermion background, while in (c) it contains all 
background contributions. 



of rejection was performed. A number of observables display separation power between 
t + t~ and qq events while being only weakly dependent on the topology classification. 
These were: the isolation angle; the event thrust; the number of neutral electromagnetic 
showers in an event; the invariant masses reconstructed using charged particles in each 
candidate r decay hemisphere of an event, and the invariant masses reconstructed using 
both charged and neutral particles in each candidate r decay hemisphere. These variables 
were used as inputs to a neural network with one hidden layer of eight neurons and a 
single output neuron. The output neuron distribution is shown in Fig. 0. A cut at 0.05 
in the output neuron gave a rejection factor of about three, leaving a background of 
0.29% from e + e~ — > qq, while rejecting only 1.3 x 10~ 3 of the remaining t + t~ events. 



Discrepancies between data and simulation in this variable are discussed in Section |6.2.1 



Background levels were estimated as accurately as possible using the data. Variables 
sensitive to a particular type of background were chosen. The background levels were 
extracted by fitting for the background contribution in these variables assuming the 
background shape from the simulation. Typically a relative precision of order 10% or 
better was achieved, ensuring systematic uncertainties on the topological branching ratios 
below the expected statistical precision. 

The background from e + e~ — > qq was estimated using the data by performing a fit 
to the distribution of the neural network output neuron. The shape of the background 
distribution was taken from simulation and the fit was performed for the normalisation 
factor. A factor 0.92 ± 0.03 was obtained compared to the simulation. 

The background due to the fully leptonic four-fermion final states was estimated from 
the data by extrapolating from the observed background in the low isolation angle region 
for events with a low visible energy, as well as for events tagged by the presence of a 
muon (for e + e~ fi + or an electron (for e + e~e + e~). The relative uncertainty obtained 
on the normalisation for these final states was ±10%. Fig. |3| shows the isolation angle 
distributions for different classes of events. 
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Figure 4: Distributions used in the study of the dileptonic background to the r + r~ pre- 
selection: a) p ra( i distribution for events containing two identified muons; b) E ra d dis- 
tribution for events containing two identified electrons. In both cases the dots are the 
data and the line is simulation. In (a) the shaded area is the background contribution 
due to the reaction e + e~ — > while in (b) it is the background contribution from 

e + e~ — ► e + e~(7) events. 



The background from dimuon events in the r + r~ sample was estimated by studying the 
momentum distribution for high momentum identified muons, and by studying the p Ta d 
distributions for events with two identified muons. This yielded a correction factor relative 
to the simulation of 0.96 ± 0.03. The p ra d distribution, after all other t + t~ preselection 
cuts, of events containing two identified muons is shown in Fig. |4]a. Similar studies 
were carried out for Bhabha background using the associated electromagnetic energy and 
anti-electron tagging in the hadron calorimeter and muon chambers. This, combined with 
studies |I7| comparing the BHWIDE and BABAMC Bhabha generators indicated that a 
rescaling upwards byl.l5±0.15of the Bhabha background estimated from simulation 
was necessary. In addition, further cross-checks were performed by studying the p ra d 
and E ra d distributions in the region outside the one occupied by r + r~ events. The E ra d 
distribution, after all other t + t~ preselection cuts, for events with two identified electrons 
is shown in Fig. |]b. 

The cosmic ray contamination was estimated from the data to be (50 ± 8) x 10 -5 by 
extrapolating from the number of events seen outside the cuts in impact parameter. 

The total background was estimated to be (1.51 ± 0.10)%, with the following break- 
down: (0.11 ± 0.01)% from fin al states; (0.40 ± 0.07)% from Bhabha events; 
(0.29 ±0.01)% from qq events; (0.27 ±0.03)% from e+e'e+e" final states; (0.10 ±0.01)% 
from e + e-/i+/i" events; (0.27 ± 0.03)% from e+e'r+r- final states; (0.02 ± 0.01)% from 
e + e~qq final states; (0.05 ± 0.01)% from cosmic rays. The background from ji + pT fi + yT , 
fi + {i~T + r~ and t + t~t + t~ final states was negligible. 

Satisfactory agreement with the simulation was observed for each year separately, and 
estimated efficiencies and backgrounds were compatible between the different years. 
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5 Topology reconstruction 

5.1 Charged particle tracking 

The charged particle track pattern recognition was based on a robust algorithm de- 
signed to minimise the number of bad associations between hits produced by different 
particles in the different elements of the tracking system. In general, most tracks contain 
hits in all the tracking subdetectors, the VD, ID, TPC and OD. Most importantly, for 
the rejection of photon conversion products, the attachment of VD hits to a track must 
be as efficient and as unambiguous as possible. The VD efficiency and misassociation 
probability in simulation have been tuned to the data with e + e~ — ► events and 



with large samples of hadronic events for the purposes of heavy quark tagging pC] . 

The TPC was the most important subdetector for charged particle reconstruction. The 
section of a track reconstructed in the TPC had to contain at least three space points out 
of a maximum possible 16. In dense topologies, where tracks lay closer together than the 
inherent two-track resolution of the detector, the clustering algorithm used to reconstruct 
space points used the information contained in the profile of a single cluster to resolve 
two tracks down to distances of typically 2 to 3 mm in r-0. Where tracks were closer 
than this, a space point was associated simultaneously to each of the tracks with which 
it was compatible and assigned a large error so as not to introduce systematic biases in 
the reconstructed parameters for these tracks. For regions away from the TPC sector 
boundaries, the probability to include either the TPC or the OD on an isolated track was 
close to 100%. 

After the data had been passed through the standard algorithm, a second algorithm 
was run which produced VD-ID tracks from combinations of hits in the VD and the ID 
jet chamber. This algorithm recuperated charged particle tracks around the boundary 
regions of the TPC or cases where the two-track resolution of the TPC prevented it from 
resolving both charged particles. The two-track resolution of the VD is more than one 
order of magnitude better in r-0 than the TPC. However the VD-ID tracks had poorly 
reconstructed momenta and particle charge sign determination due to their short length. 
This had no consequences for the charged particle multiplicity reconstruction. 



5.2 Photon conversions 

A good understanding of the photons converting in material before the tracking sub- 
detectors is important in order to measure correctly the true charged particle multiplicity 
in a t decay. About 7% of photons interact with the material before the TPC gas vol- 
ume, giving an e + e~ pair detected in some of the tracking chambers. In particular, 
unreconstructed conversions form an important potential background from lower multi- 
plicity topologies in higher multiplicity samples. Furthermore, to measure in data both 
the conversion reconstruction efficiency and the amount of material, it was necessary to 
study the rate of both reconstructed and unreconstructed converted photons. 

Converted photons were reconstructed using the algorithms described below. These 
are described in more detail in PT1. Photons converting in the beam-pipe, VD, ID and 



TPC inner wall were reconstructed using charged particle tracks observed in the TPC. 
For each TPC track the position of the tangent to the helix which passed through the in- 
teraction region was calculated. This point gave the estimated position of the conversion. 
All tracks for which this point was compatible with the interaction region within one 
standard deviation were neglected as conversion candidates. If two conversion candidates 
were found with compatible conversion points, they were accepted as a converted photon 
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provided that: the two tracks had opposite curvature; the mean conversion radius of the 
two tracks was greater than 5 cm; at least one of the tracks had no associated point 
inside the mean conversion radius; the angles of the two tracks at the conversion radii 
agreed within 30 mrad in <fi and 15 mrad in 9. This algorithm had an efficiency, esti- 
mated from simulation, of about 60% for photon energies below 6 GeV falling to about 
20% at 10 GeV. It was complemented by an algorithm using a fit to a joint origin for 
two charged particle tracks assuming zero opening angle, with similar requirements on 
the position of the reconstructed conversion radius to the previous algorithm. This had 
an efficiency of about 40% for photon energies above 8 GeV. The average efficiency for 
the combination of both algorithms was (68.1 ± 0.2)% for the photon energy spectrum in 
simulated 1-prong r decays and (59.8 ± 0.4)% for 3-prong decays. 

In 3-prong decays, the probability of reconstructing a it + it~ pair incorrectly as a 
photon conversion was estimated from simulation to be 1 x 10~ 4 . For systematic studies 
the uncertainty on this was taken as 100% of the size of the effect. The probability of 
reconstructing a Kg — > 7r + 7r~ decay incorrectly as a photon conversion was estimated 
from simulation to be 1.5%, however this did not alter the observed charged particle 
multiplicity. 

The rates of reconstructed and unreconstructed converted photons were studied as 
a function of the radial coordinate of the estimated conversion point. The conversions 
were separated into those taking place in the beam-pipe, VD and ID jet chamber inner 
wall, the ID jet chamber sensitive volume, the ID jet chamber outer wall and ID Trig- 
ger Layers, and the TPC inner wall. The bulk of the material was in the TPC inner 
wall. In order to be as independent as possible of the t + t~ presample and modelling of 
the photon production rate in r decays, the study was performed using e + e~ — > /i + (j,~~f 
and e + e~ — > e + e~7 radiative events. The events were selected using electron or muon 
identification and used kinematic constraints to predict the direction and energy of the 
photon. The rates were recorded for events with an HPC shower, reconstructed conver- 
sion, or no reconstructed conversion but some charged particle tracks consistent with an 
electron-positron pair arising from a converted photon. In the latter case the radius of 
the conversion point was taken as being that of the material layer lying just inside the 
innermost first measured point on either of the two tracks believed to be due to the con- 
version. The rates were normalised to the total number of observed photons, including 
conversions. The identical procedure was applied to the simulation. Correction factors 
for the reconstructed conversion rate were obtained from the ratio A c °" v of the rates R7 C 

rec l* 

and R r s ec of reconstructed conversions in data and simulation respectively. Analogous 
correction factors for the unreconstructed conversion rate were obtained from the ratio 

. conv r | j-jUnrec , j-^unrec r -. 

A of the rates H r , and K„ of unreconstructed conversions. 

unrec <* * 

The measured correction factors were used to estimate the reconstruction efficiency 
and amount of material. The ratio A c ° nv is given by 

rec ^ ' J 

,-jrec j-jQ ^EM conv 

conv _ K d _ ^ X d ^ . . 

reC R reC J^P J^EM ^conv ' V / 

The subscripts d and s stand for data and simulation respectively. R d s is the total rate of 
photons, X d ™ is the probability of a conversion taking place (proportional to the amount 
of material in terms of radiation lengths) and e d °™ is the reconstruction efficiency. The 

. conv . . , 

ratio A is given by 

unrec - J 

jyunrec „q ^EM , conv 

conv _ K d _ ^ X d ^ 1 — e d , . 

unrec j^unrec j^q j^EM ^ ^conv * V U / 



11 



Radius 


Conversions 


Lower 


Upper 


. conv 
rec 


. conv 
unrec 


conv i conv 


X d l X s 


Ocm 


13.0cm 


1.00±0.03 


1.14±0.05 


0.96±0.03 


1.04±0.03 


13.0cm 


22.0cm 


1.34±0.05 


1.70±0.12 


0.93±0.06 


1.45±0.05 


22.0cm 


28.5cm 


1.23±0.03 


1.50±0.09 


0.94±0.05 
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Table 1: Factors used to correct simulated photon conversions and data/ simulation dis- 
crepancies as a function of the radius of the interaction point. Column 5 shows the 
data/simulation ratio of reconstruction efficiency and column 6 the equivalent ratio for 
the material. 



The conversion reconstruction efficiency in data e° d ° can be obtained from the data alone 
via the relation 

rec 

conv / t\ 

7~> rec I jjunrec •> \ ) 

K d + K d 

independently of any assumptions on the material or, if Rf d c and R d " rec are measured 
on the same data sample, on the incident photon flux. If the ratio R d / R° s is known, 

P , conv , conv , conv , . -. . , conv r , . , 

measurements of A , A and e w , combined with e c from the simulation, can 

rec ' unrec u > a > 

E Aif F! M~ 

provide an estimation of the material ratio X d /X s for electromagnetic interactions. 
The ratio R d /R° s of the rates of photons produced in the e + e~ — > fi + fi~'y and e + e _ — > 
e + e~7 radiative events was taken to be unity with a systematic uncertainty of ±2% |22|. 



The factors A and A obtained as a function of interaction radius are given in 

rec unrec ^ J 

Table [3]. These factors were used to reweight the simulated events used in the estimation 
of the selection efficiencies to account for data/simulation discrepancies in the material 
and reconstruction efficiency. The distributions of the conversion radius and energy for 
reconstructed converted photons are shown in Fig. [|, where the simulation has been cor- 
rected with the factors obtained in the e + e~ — > / u + / u _ 7 and e + e~ — > e + e~7 events. There 
are still some localised discrepancies between data and simulation due to simplifications 
in the detector material description in the simulation program. For the r presample, the 
rate of reconstructed conversions observed in data was 1.004±0.008 times that predicted 
by the simulation after the application of the correction factors, in excellent agreement. 

The data taken in different years were studied separately, in particular before and 
after the beginning of 1995, the point after which the material in the ID trigger layers 
was substantially reduced. Consistency was observed between the different years and the 
change in conversion rate observed in the ID trigger layers between 1992-1994 and 1995 
was consistent with the expectation. 

The conclusion, integrating up to a radius of 50 cm, is that the material is underesti- 
mated in the detector simulation by about 10% in terms of radiation lengths, while the 
conversion reconstruction algorithm efficiency is similar in data and simulation. 

The kinematical conversion reconstruction was complemented by conversion rejection 
based on an electron identification algorithm which combined information from the 
TPC ionisation and the HPC electromagnetic shower reconstruction. This algorithm 
also rejected many tt° Dalitz decays, which was particularly important for rejecting back- 
grounds in the 5-prong r decay sample. A conversion was flagged if there was at least 
one tightly tagged electron in a hemisphere with greater than one track in the TPC or 
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Figure 5: a) Distribution in e + e~ — > t + t - events of the estimated conversion radius for 
photons converted in front of the TPC. Data are dots, the line is simulation after applying 
correction factors from e + e _ — > /i + |U~7 and e + e~ — > e + e~7 events, b) Reconstructed 
conversion photon energy in e + e~ — > t + t~ events. Data are dots, the line is simulation 
after applying correction factors. 



if there were a pair of oppositely charged particles both of which were loosely tagged 
as electrons. These criteria were chosen to minimise the probability of misidentifying a 
7r + 7r _ pair as e + e~ rather than to maximise the conversion identification efficiency In 
simulation, the efficiency to reject the remaining conversions was (59.6±0.6)% in 1-prong 
decays and (54.5±0.8)% in 3-prong decays. The electron-positron pair rejection also 
identified (72.7±1.2)% of Dalitz pairs from tt° decays in simulated 1-prong decays, giving 
an improved classification for these decays. The efficiency was (68.3±3.1)% for 3-prong 
decays. The efficiency to identify Dalitz pairs was higher than for conversions due to 
biases arising from the kinematic conversion reconstruction algorithm which had already 
identified and reconstructed most converted photons. 

The conversion rejection efficiency of the electron identification algorithm was cross- 
checked in t decays containing three charged particle tracks, from the rate of oppositely 
charged pairs of particles which satisfied the electron tagging requirements outlined above 
and where neither particle had a hit in the VD. The ratio of the efficiency in data to that 
in simulation was estimated to be 0.96±0.05, in good agreement with unity. The relative 
systematic uncertainty on the efficiency of this algorithm was conservatively taken to 
be 10%, to account for any systematic effects in the estimation of the ratio and in its 
application to all classes of charged particle tracks. 

Fig. |6| shows the invariant mass distribution for pairs of oppositely charged particles 
in r decay hemispheres containing three charged particles, where one or both of the pair 
were tagged as an electron candidate and where both tracks had hits in the VD. The 
electron mass was assumed for both particles. The VD association requirement rejected 
most photon conversions, which occurred after the VD, while not affecting the signal of 
7T° Dalitz decays, or of hadrons misidentified as electrons. The contribution from Dalitz 
decays peaks at low masses, below the 7r° mass. The tail from Dalitz decays in the 
region of masses around 200 MeV to 600 MeV arises from decays where only one of the 
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Figure 6: Invariant mass m e + e - distributions for pairs of oppositely charged particles 
with one or both particles identified as electron candidates and where both tracks have 
hits in the VD. Data are dots, the simulation is the line. The shaded area shows the 
contribution from 7r° — > 7e + e~ decays. The region m e + e - > 200 MeV is dominated by 
misidentified -n + n~ pairs. 



electron or positron was identified and there were two possible combinations of oppositely 
charged particles. In such decays only half the possible combinations were due to the e + e~ 
pair. In these cases each combination has been given a weight of one half. The region 
of mass greater than 200 MeV was dominated by 7r + 7r~ pairs where one of the pions 
was misidentified as an electron. This permitted a determination from the data of the 
probability to misidentify a 3-prong r decay as a 1-prong decay containing a conversion 
candidate. The misidentification probability was estimated to be (1.29 ± 0.09)%. The 
uncertainty was dominated by data statistics and contained a contribution from the VD 
association efficiency From the ratio of data to simulation in the low mass spike it was 
estimated that the efficiency to reject a tt° Dalitz decay in data was 0.90±0.07 times 
the efficiency obtained in simulation. The uncertainty included contributions from the 
error on the Particle Data Group best estimation of the 7r° — > 7e + e~ branching ratio ||, 
the uncertainty in the VD efficiency and the limited simulation statistics in addition to 
the statistical uncertainty. No correction was applied to the simulation and a relative 
uncertainty of 10% on the Dalitz rejection efficiency was taken. 

Most 5-rays were rejected by the requirement of VD hits on a track. Events in simula- 
tion containing a 5-ray were reweighted with the radius dependent factors describing the 
ratio of material in terms of radiation lengths X d / X s in Table |I|. There is some am- 
biguity in the choice of this factor. The associated systematic uncertainties are discussed 
in Section |6.2.3 . 

5.3 Nuclear reinteractions 

About 3% of hadrons reinteract inelastically with the material before the TPC gas 
volume, and produce typically up to 10 charged secondary particles. These reinteractions 
were reconstructed by an algorithm which was designed to find secondary reinteraction 
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vertices using the tracks from outgoing charged particles produced in nuclear reinterac- 
tions. 

The algorithm was robust with respect to the number of charged particles produced 
and to the position of the reinteraction with the DELPHI detector material. Only tracks 
which had impact parameters inconsistent with coming from the interaction region were 
considered. The algorithm iterated over all pairs of tracks, and attempted to fit a ver- 
tex in three dimensions if neither of the tracks had already been included in an earlier 
successfully fitted vertex. However, if one of the tracks had already been included in 
an earlier fitted vertex, an attempt was made to include the other track in that vertex. 
This was the main part of the algorithm permitting the reconstruction of vertices with an 
arbitrary number of outgoing charged particle tracks. In the rare case where both tracks 
had already been associated to different fitted vertices, an attempt was made to merge 
those two vertices, if their positions were compatible. In all candidate reinter- 

action vertex had to pass a x 2 cu ^ an d all outgoing tracks had to have a first observed 
hit consistent with the particle originating from the reconstructed vertex. Candidate 
vertices with two outgoing tracks and kinematically compatible with the decay of a K® 
or A produced in the interaction region were not considered. 

A search was then performed for the incoming particle which caused the interaction. 
Tracks which had been reconstructed in only the VD or the VD and ID were extrapolated 
to the candidate nuclear reinteraction vertex, and associated to the vertex if consistency 
was found. This prevented double counting of VD-ID tracks and reconstructed nuclear 
reinteractions. If there was no such track the momentum and charge of the incoming 
hadron were estimated from the sum over the outgoing particles associated to the vertex. 
For t decays it is a good approximation to assume all such nuclear reinteraction vertices 
are caused by charged particles, as K° mesons are the only source of neutral hadrons 
producing nuclear reinteractions and are infrequently produced. 

In simulation, the algorithm was (82.4±0.2)% efficient in 1-prong r decays for incoming 
particle momenta above 1 GeV/c, and (83.2±0.4)% efficient in 3-prong r decays. The 
efficiency had little dependence on the r decay multiplicity or on the charged secondary 
particle multiplicity in the nuclear reinteraction. The position of the reinteraction vertex 
was reconstructed with a precision of better than 1 mm in three dimensions. 

In a similar manner as for the converted photons, the reconstructed and unrecon- 
structed nuclear reinteractions were studied as a function of the reinteraction radius. 
The rates of reconstructed nuclear reinteractions were obtained from the r + r~ sample 
in data and simulation. The radial distribution of reconstructed hadronic reinteractions 
is shown in Fig. [7|. Some localised discrepancies, more marked than in the case of the 
photon conversions, are visible in the detailed description of the ID trigger layers and 
TPC in the simulation. The flat tail for radii greater than 35 cm is due to reinteractions 
in the TPC cathode plane at 9 ~ 90°. 

The unreconstructed nuclear reinteractions were measured by studying the multiplicity 
distribution for charged particle tracks with impact parameters in the r-0 plane greater 
than 1.5 cm in the hemisphere opposite a well reconstructed leptonic r decay candidate. 
Simulation studies showed that this was dominated by r decays with an unreconstructed 
nuclear reinteraction for multiplicities of four and greater. The r decay hemispheres used 
to estimate this effect had to have no reconstructed nuclear reinteraction. This removed 
the risk of a false measurement arising from potential discrepancies in the charged particle 
multiplicity due to data/simulation differences in the track inclusion efficiency of the 
nuclear reinteraction reconstruction algorithm. The reinteraction radius was taken as the 
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Figure 7: Distribution in e + e~ — ► t + t~ events of reconstructed radius of the inelastic 
nuclear reinteractions. Data are dots, the line is simulation after applying correction 
factors. 

radius of the material surface lying immediately inside the first measured point on any 
of the charged particle tracks with a high impact parameter. 

The measured correction factors A nucl and A nucl , analogous to the quantities A c °" v and 

rec unrec ' 1 rec 

A° for conversions, are given in Table 0, together with the data/simulation ratio of the 
efficiencies for reconstruction of inelastic nuclear reinteractions / e"" c and the derived 
data/simulation ratio of material in terms of nuclear interaction lengths X d ° /X™ uc . The 
region of the ID jet chamber gas volume, (13 cm to 22 cm) has been included together with 
the region at lower radius as it has a very small number of nuclear interaction lengths, of 
order 10~ 4 . These results indicate that for nuclear interactions the material at lower radii 
than the TPC gas volume is over-estimated by about 9% in the simulation program while 
the reconstruction efficiency is consistent in data and simulation. There is a significant 
difference in the estimated data/simulation material ratios for nuclear and electromag- 
netic interactions. This can be attributed to the complexity of the composite material 
structure together with the means by which each of these quantities is introduced in the 
simulation program, as well as to possible weaknesses in the modelling of the hadronic 
reinteractions in the simulation program. The same effects are observed in e + e~ — > qq 
events. The uncertainty in the initial production rate of hadrons in r decays is relatively 
small and does not give an important systematic uncertainty in the measurement of A uc 
and A 

unrec 

As for the photon conversions, the nuclear reinteractions were studied for each year 
separately, particularly before and after the upgrade of the ID trigger layer in 1995. 

Simulation studies indicate that the probability of reconstructing a fake nuclear rein- 
teraction from r decay secondary particles where none was present was 5 x 10 -4 in 3-prong 
t decays. However most such events were still classified as 3-prongs. The associated un- 
certainty was taken to be 100% of the size of the effect. The probability of reconstructing 
a fake nuclear reinteraction from a K§ decay was negligible. 

The probability that a hadron produced in a r decay would undergo an elastic nu- 
clear reinteraction before the TPC sensitive volume was estimated from simulation to 
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Inelastic nuclear reinteractions 



Lower 


Upper 


. nucl 
rec 


. nucl 
unrec 


nucL 1 nucl 


-r^nucl i -tunnel 


Ocm 


22.0cm 


0.91±0.08 


1.01±0.02 


0.95±0.05 


0.96±0.04 


22.0cm 


28.5cm 


0.93±0.04 


1.36±0.15 


0.96±0.13 


0.97±0.04 


28.5cm 


50cm 


0.86±0.03 


1.00±0.03 


0.98±0.03 


0.87±0.03 


Ocm 


50cm 


0.89±0.02 


1.01±0.02 


0.98±0.02 


0.91±0.02 



Table 2: Factors used to correct simulated nuclear reinteractions for data/ simulation 
discrepancies as a function of the radius of the reinteraction point. Column 5 shows the 
data/simulation ratio of reconstruction efficiency and column 6 the equivalent ratio for 
the material. 



be (3.29±0.05)%. The typical signature of such an interaction was a charged particle 
track with a kink at the scattering point in the detector material. For large scatters the 
reconstruction algorithm would not reconstuct one but two tracks, one before and one 
after the scatter. This led typically to a signature of one VD-ID track and one TPC-OD 
track, and did not alter the observed track multiplicity in the TPC and OD. 



5.4 Topology selection criteria 

By judicious application of the VD hits requirement, and taking into account that the 
r must decay into an odd number of charged particles, algorithms were constructed to 
minimise the sensitivity to the VD association efficiency and the knowledge of the mate- 
rial reinteractions; these are described below. In the following a "good" track is defined 
as a charged particle track which has associated TPC or OD hits and is not identified 
as originating from a conversion or nuclear reinteraction in the detector material. For 
simplicity, the VD-ID track classification, in addition to including tracks reconstructed 
using only the VD and ID detectors, includes tracks reconstructed as the ingoing charged 
hadron to a nuclear reinteraction. In r decay hemispheres where the presence of a con- 
verted photon was signalled by the electron identification algorithm, the observed number 
of good tracks was reduced by two (or by one in the rare case of an even number of good 
tracks) for the purposes of the topology classification. 

Individual r decay hemispheres were allotted to five different classes dependent on the 
number of reconstructed tracks in the TPC, VD, etc. The first three of these classes 
correspond closely to events with clean reconstruction of a 1-prong, 3-prong or 5-prong 
topology. The other two classes contain the small number of r decay hemispheres where 
there was more ambiguity in the reconstruction of the topology, but where some discrim- 
ination was still possible. 

A 1-prong r decay was defined as a r decay hemisphere satisfying any of the following 
criteria: 

• only one good track with associated VD hits, and no other tracks with associated 
VD hits; 

• only one good track, without VD or ID hits, and one VD-ID track; 

• no good tracks, and only one VD-ID track. 

3-prong r decays were isolated by demanding r decay hemispheres satisfying at least one 
of the following sets of criteria: 
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• three, four or five good tracks, of which either two or three had VD hits; 

• two good tracks with associated VD hits, plus one VD-ID track; 

• one good track with associated VD hits, plus one or two VD-ID tracks pointing 
within 3° in azimuth of a TPC sector boundary. 

Candidate 5-prong r decays were selected if they satisfied at least one of the following 
topological criteria: 

• five good tracks of which at least four had two or more associated VD hits; 

• four good tracks with VD hits, and one other VD-ID track. 

Additional criteria were applied in the selection of 5-prong r decays due to the large 
potential background from hadronic Z decays and mis-reconstructed 3-prong r decays. 
The background originating from v^h^ > In final states with a Dalitz decay was 
expected to occur at a similar level to the signal. The electron rejection criteria described 
in Section |0] reduced this background by about 70%, and it was further suppressed by 
requiring that all good tracks had a momentum greater than 1 GeV/c. To reject Z —>■ qq 
events it was required that the total momentum of the 5-prong system be greater than 
20 GeV/c. Only good tracks were included in the calculation of this quantity. 

These three classes accounted for 97.6% of candidate r decays in the r + r~ sample. 
The remaining 2.4% of candidate r decays were mostly 1-prong and 3-prong r decays 
with some pattern recognition failure or detector inefficiency. These were classified into 
two categories, 1' and 3', corresponding to 1-prong and 3-prong r decays respectively. 
Candidate r decay hemispheres of type 1' had to satisfy at least one of the following 
criteria: 

• one good track, with associated VD hits, plus two VD-ID tracks not pointing within 
3° in azimuth of a TPC sector boundary; 

• one good track, plus a conversion candidate selected with the electron identification 
algorithm. 

The 3' category contained all other r decay hemispheres failing to pass any of the previous 
classification requirements. In particular it contained those events with either four or five 
good tracks of which either four or more had associated VD hits, but not satisfying the 
5-prong class requirements. This criteria had a higher efficiency for 5-prong decays than 
for 3-prong decays, but due to the much greater 3-prong branching ratio the 3' category 
was dominated by 3-prong decays. Table ^ contains the efficiencies of these selection 
requirements for the different exclusive r decay modes and the inclusive single-hemisphere 
topological selections, as obtained from simulation and after corrections for observed 
discrepancies between data and simulation in the rate and reconstruction efficiency of 
material reinteractions. 

Events were classified according to both the hemisphere classifications into 14 event 
topology classes: 1-1, 1-1', 1-3, 1-3', 1-5, l'-l', l'-3, l'-3', l'-5, 3-3, 3-3', 3-5, 3'-3' and 3'-5. 
The class 5-5 was not included as these events were removed by the multiplicity cut in 
the t + t~ preselection and in any case would give a negligible signal contribution. The 
selection of event classes up to multiplicity five assumes that the inclusive branching ratio 
of the r to seven charged particles is negligible. While inclusion of higher multiplicity 
classes is possible, the DELPHI sample size is insufficient to reach the level of precision 
obtained by CLEO [T§. 



The various sources of non-r+r background are detailed in Table |] for each of the 
event topology classes. 
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true t 




Topology Classification 


decay mode 


selection 


1 


3 


5 


V 


3' 




50.60±0.07 


99.95±0.00 


0.00±0.00 


0.00±0.00 


0.02±0.00 


0.02±0.00 


\TVrVy, 


53.31±0.07 


99.96±0.00 


0.00±0.00 


0.00±0.00 


0.01±0.00 


0.03±0.00 


7T~V T 


49.69±0.09 


99.88±0.01 


0.04±0.01 


0.00±0.00 


0.00±0.00 


0.08±0.01 


7T~7T°l> T 


51.77±0.06 


97.87±0.03 


0.60±0.01 


0.00±0.00 


0.84±0.02 


0.69±0.01 




51.07±0.11 


95.88±0.06 


1.25±0.03 


0.00±0.00 


1.42±0.04 


1.45±0.04 


7T~37T°b> T 


48.89±0.25 


94.36±0.16 


1.68±0.09 


0.00±0.00 


1.83±0.10 


2.13±0.10 


K~v T 


49.43±0.36 


99.90±0.03 


0.02±0.02 


0.00±0.00 


0.01±0.01 


0.06±0.03 


K-n°u T 


51.40±0.47 


97.66±0.20 


0.85±0.12 


0.00±0.00 


0.80±0.12 


0.69±0.11 


K-27i°u T 


50.42±1.12 


94.65±0.71 


2.28±0.47 


0.00±0.00 


1.39±0.37 


1.68±0.40 




53.10±0.48 


99.79±0.06 


0.07±0.03 


0.00±0.00 


0.00±0.00 


0.14±0.05 


7r-^ L 7rV 


51.85±0.73 


97.32±0.33 


0.78±0.18 


0.00±0.00 


0.78±0.18 


1.11±0.21 




54.60±0.87 


99.78±0.11 


0.11±0.08 


0.00±0.00 


0.00±0.00 


0.11±0.08 




52.66±1.24 


96.71 ±0.61 


0.94±0.33 


0.00±0.00 


0.94±0.33 


1.41±0.40 




52.82±1.04 


95.12±0.62 


3.72±0.54 


0.00±0.00 


0.08±0.08 


1.07±0.30 




52.17±0.48 


94.48±0.30 


4.30±0.27 


0.00±0.00 


0.16±0.05 


1.06±0.14 




50.78±0.73 


92.64±0.54 


4.65±0.43 


0.00±0.00 


0.55±0.15 


2.16±0.30 


K-K s °u T 


52.38±0.86 


94.50±0.54 


4.42±0.49 


0.00±0.00 


0.06±0.06 


1.02±0.24 


K-Kg7r°u T 


51.32±1.32 


92.56±0.97 


5.01±0.80 


0.00±0.00 


0.41±0.23 


2.03±0.52 


n-2K°v T 


46.34±1.80 


86.72±1.80 


10.45±1.63 


0.00±0.00 


0.00±0.00 


2.82±0.88 


1-prong 


51.42±0.03 


98.71±0.01 


0.42±0.01 


0.00±0.00 


0.44±0.01 


0.43±0.01 


27T~7T + U T 


54.71±0.11 


0.90±0.03 


90.26±0.09 


0.01±0.00 


2.10±0.04 


6.74±0.07 


27T _ 7r + 7T Z/ T 


53.88±0.13 


1.26±0.04 


86.39±0.12 


0.10±0.01 


3.15±0.06 


9.10±0.10 


37r ± 27r°z/ r 


53.14±0.46 


1.37±0.15 


83.64±0.46 


0.22±0.06 


3.68±0.24 


11.09±0.39 


37r ± 37r°z/ r 


52.13±1.06 


1.46±0.35 


78.73±1.20 


0.17±0.12 


4.74±0.62 


14.90±1.05 


K~7T~7T + V T 


54.64±0.56 


1.03±0.15 


90.35±0.45 


0.00±0.00 


1.70±0.20 


6.92±0.38 


K-K + ti + v t 


53.87±0.90 


2.08±0.35 


87.23±0.82 


0.00±0.00 


1.10±0.26 


9.59±0.73 


3-prong 


54.32±0.08 


1.06±0.02 


88.51±0.07 


0.05±0.00 


2.54±0.03 


7.84±0.06 


?>'K~2-K + V T 


49.63±1.19 


0.11±0.11 


12.63±1.13 


57.52±1.67 


0.23±0.16 


29.51±1.55 


57r ± 7r°z/ T 


48.91±2.23 


0.00±0.00 


15.04±2.28 


52.85±3.18 


0.81±0.57 


31.30±2.96 


5-prong 


49.47±1.05 


0.09±0.09 


13.16±1.01 


56.49±1.48 


0.36±0.18 


29.90±1.37 



Table 3: Estimates of the selection and topology classification efficiencies, in percent, 
for different exclusive decay modes, as obtained from simulation. The efficiencies are 
corrected for observed discrepancies between data and simulation in the rate and recon- 
struction efficiency of material reinteractions. The quoted uncertainties are from the 
simulation statistics only. 

6 The fit and systematics 

6.1 Fitting procedure 

A reweighting technique was used to take into account the observed data/simulation 
discrepancies in the rates and reconstruction efficiencies of nuclear reinteractions, photon 
conversions and 5-rays. Each secondary or tertiary particle i in a simulated r + r~ event k 
was given a weight W[ art . For particles interacting with the detector material this weight 
was obtained from the studies described in Section |5] and summarised in Tables [1] and |2|. 
These weights were a function of the radius and type of the reinteraction, as follows: 
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Source of 


All 


Event Topology 


Background 


Topologies 


1-1 


1-1' 


1-3 


1-3' 


1-5 


l'-l' 


l'-3 


T-3' 


l'-5 


3-3 


3-3' 


3-5 


3'-3' 


3'-5 




0.11±0.01 


0.19 


0.06 


0.00 


0.07 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


e + e~ 


0.40±0.07 


0.65 


0.00 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


qq 


0.29±0.01 


0.03 


0.00 


0.32 


3.39 


1.01 


0.00 


1.16 


0.00 


0.00 


1.61 


14.8 


21.2 


77.0 


0.00 


e + e~e + e~ 


0.27±0.03 


0.49 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


e + e~ fi + fi~ 


0.10±0.01 


0.16 


0.00 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.33 


0.00 


0.00 


0.00 


e + e~T + T~ 


0.27±0.03 


0.35 


3.91 


0.23 


0.76 


0.00 


21.1 


0.97 


2.88 


0.00 


0.20 


0.60 


0.00 


0.45 


0.00 


e + e~qq 


0.02±0.01 


0.03 


0.00 


0.03 


0.44 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


cosmic rays 


0.05±0.01 


0.07 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


Total 


1.51±0.10 


1.97 


3.98 


0.60 


4.66 


1.01 


21.1 


2.13 


2.88 


0.00 


1.81 


15.7 


21.2 


77.4 


0.00 



Table 4: Fractional non-r + r~ backgrounds (in %) in the t + t~ sample and the different 
event topology classes. The classes l'-5 and 3' -5 have very small populations and the 
estimated background is zero due to statistical fluctuations in the simulation. 



• reconstructed converted photons: Wf art = A coni '(r), 

• unreconstructed converted photons: Wf art = A c ° nv (r) 

• reconstructed inelastic nuclear reinteractions: Wf art = A" uc (r) 

' rec ^ ' 

• unreconstructed inelastic nuclear reinteractions: Wf art = A"" c (r); 

• 5-rays: ^r rt = Xr7< M W; 

• bremsstrahlung emitting electrons: W[ art = X d /X s (r); 

• elastic nuclear reinteractions: Wf art = x™ 01 / X™ ucl (r) . 

It was also necessary to reweight particles which did not reinteract so as to maintain 
the internal normalisation of the simulation sample with respect to parameters such 
as angular distributions, momentum distributions, or the r exclusive branching ratios. 
Failure to achieve this could lead to biases within the simulation sample as a function 
of photon multiplicities, charged hadron multiplicities, angular regions with different 
amounts of material, or could bias the preselection of the r + r~ sample with respect to 
certain types of decay mode. Non-interacting photons were given a weight Wf art given by 



w part 



n " - e w r rt 

jgint. photons 



where N 1 is the total number of photons produced by r decays or by interactions with 
material of r decay products in the simulation sample. N? nt is the number of photons 
with a material reinteraction. The sum is over all photons with a material reinteraction. 
To maintain the angular and momentum distributions the above weight calculation was 
performed separately for individual bins in a three-dimensional space of the photon polar 
angle, azimuthal angle and momentum. A similar procedure was used to reweight non- 
interacting hadrons, particles not emitting a 5-ray (a very small correction) and electrons 
without any bremsstrahlung emission. 

With every secondary and tertiary particle i in a simulated event k given a weight 
Wf art , the weight W^ vent of the event was given by the product of all these weights: 

W event = JJ w part _ (g) 
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Class 


Observed 


Fit Output 


x 1 


1 -1 


56219 


56149.0 


0.1 


1 -1' 


858 


871.0 


0.2 


1 -3 


18681 


18813.5 


0.9 


1 -3' 


2350 


2331.5 


0.1 


1 -5 


94 


95.6 


0.0 


l'-l' 


4 


3.9 


0.0 


l'-3 


131 


134.8 


0.1 


l'-3' 


16 


14.7 


0.1 


l'-5 





1.2 


1.2 


3 -3 


1481 


1451.5 


0.6 


3 -3' 


409 


357.4 


7.4 


3 -5 


17 


13.8 


0.7 


3'-3' 


76 


97.7 


4.8 


3'-5 


1 


1.5 


0.2 


All 


80337 


80337.1 


16.4 



Table 5: The second and third columns contain the number of observed events in each 
class and number predicted from the maximum likelihood fit of the r topological branching 
ratios. The fourth column contains the \ 2 contribution for each class to the overall \ 2 . 

The numbers of selected events in each event topology class are shown in Table |5]. 
A maximum likelihood fit assuming Poissonian probabilities was performed to the 
reweighted numbers of events estimated using Eqn. [I| The number of r + r~ events (N TT 
in Eqn. |3|) and the branching ratios B\ and B 5 were allowed to vary in the fit while B 3 
was constrained by the relation B\ + B 3 + B 5 = 1. The output of the fit is also shown in 
Table |. The results of this fit were: B x = (85.316 ± 0.093)%; B 3 = (14.569 ± 0.093)%; 
-B5 = (0.115±0.013)%, where only the statistical error from the fit is quoted. An estimate 
of the consistency of the fit was made by calculating a \ 2 - This took into account only 
the statistical uncertainties. It gave a x 2 / n -d.f. of 16.4/11, indicating good consistency. 
The contribution to the \ 2 horn each class is shown in Table [5| 

6.2 Systematics 

In general, the systematic uncertainties on the topological branching ratios due to 
any particular effect were estimated simultaneously for B\, B3 and B§ by repeating the 
analysis, including the t + t~ preselection, after modifying the relevant variable in the 
simulation. This accounted for correlations in the systematic uncertainties between the 
different branching ratios and correlations in efficiencies and backgrounds between the 
different event classes and the t + t~ preselection. 

6.2.1 Preselection 

The main systematic effects of the t + t~ selection criteria on the result can arise 
through the mis-calibration of the quantities used in the selection. These quantities 
can be classified into cuts related to energy or momentum measurements, such as E ra ^ 
p ra d and E vis , or those cuts related to multiplicity, such as the isolation angle, which 
can be incorrectly estimated if extra secondary particles are produced in reinteractions 
with detector material. This second effect is taken account of by the systematics in 
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the reweighting of the secondary interactions, discussed in Section |6.2.3j . A separate 
systematic uncertainty is included for the energy and momentum scales below. Other 
effects such as tracking efficiency, trigger efficiency and uncertainties on the r exclusive 
branching ratios can have an effect on the r + r~ preselection efficiency. These sources of 
systematic uncertainty are discussed below, and in general the systematic uncertainties 
on the topological branching ratios take into account effects in the preselection. 

Any remnant discrepancies due to r + r~ preselection criteria were checked by studying 
the agreement between data and simulation for the distributions of the t + t~ selection 
variable for each of the different event topology classes. With the corrections to the 
background levels discussed in Sections [4.2| and [5], good agreement was found, in particular 
in the regions of the cuts. 

The uncertainties due to the backgrounds from non-r + r~ sources were estimated by 
varying the background normalisations by their uncertainties obtained as described in 
Section |J and are listed in Table |6|. The distribution of the neural network variable used 
for qq rejection (see Fig. |2|) displays a 15% excess in the data compared with simulation 
in the region [0.05;0.8] of the output neuron distribution. This region was dominated by 
t + t~ events and was studied to see if there were any discrepancies in the input variables, 
and by checking the multiplicity distributions of the events. No discrepancies were found. 
The systematic uncertainty was estimated by rescaling the numbers of rejected signal 
events in all classes by 15%. 



6.2.2 Tracking 

Within the angular acceptance of this analysis, there are four tracking detectors which 
can contribute to the track reconstruction of all charged particles: the VD, ID, TPC and 
OD. The reconstruction of "good" tracks is strongly dependent on the TPC with its 
full three-dimensional readout. The redundancy in the track reconstruction obtained by 
the inclusion of r decay hemispheres with only a VD-ID track in the r + r~ sample and 
1-prong subsample reduces greatly the sensitivity of the measurement to inefficiencies 
within a given subdetector, and allows direct cross-checks to be made. 

In simulation, for non-interacting particles, the efficiency for the TPC to reconstruct an 
isolated charged particle which passes through a sensitive gas volume far from dead regions 
is 99.95%. However, even with the redundancy of the different tracking subdetector 
components in DELPHI, to measure this in data directly is difficult because of material 
reinteractions between different subdetectors which can cause the particle to be lost before 
entering the TPC sensitive volume. These effects tend to reduce the measured efficiency 
For the data, the efficiency of the TPC in its sensitive regions was estimated using the 
redundancy of the tracking system to be (99.35 ± 0.05)%. The identical procedure was 
applied to the simulation, yielding an efficiency of (99.31 ±0.01)%, in excellent agreement 
with the data, and implying that the modelling of the TPC efficiency in the simulation was 
accurate. The true inefficiency in simulation is a factor of 12 lower than the inefficiency 
estimated by this method. The uncertainty on the TPC reconstruction efficiency in the 
TPC sensitive region was taken conservatively to be 0.05% to account for any systematic 
effects in its estimation. 

The TPC efficiency for isolated tracks in the sector boundary region was studied using 
inclusive low multiplicity events but excluding the t + t~ events. This sample consists 
of radiative dilepton events and the high energy part of the two-photon leptonic event 
spectrum. The number of reconstructed tracks in the TPC boundary region was compared 
to that expected by normalisation of the sensitive region of the TPC sectors, in bins of 
momentum. The loss of tracks was compatible with zero below about 5 GeV/c, rising to 
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3.5% at 45 GeV/c. Data and simulation agree well in their behaviour. The inefficiency 
estimated from data for the t + t~ momentum distribution was (2.73±0.04)%, compatible 
with the result in simulated r + r _ events. The TPC efficiency was varied by throwing 
away tracks containing TPC hits in simulation and repeating the analysis, including 
the t + t~ preselection stage. The error was scaled by a factor two to account for any 
systematic effects in the procedure. 

The attachment of VD hits to a track is the main criterion used in this analysis to 
determine the multiplicity of a r decay. The association efficiency and mis-association 
probability of a VD hit to be attached to a charged particle track has been studied [2(J 
on large samples of e + e~ — > qq events, and on isolated topologies such as final 
states. These studies indicate that the efficiency in simulation and data agree to within 
±2%. The associated systematic uncertainties were obtained from the shifts in the results 
observed when randomly removing 2% of associated VD hits and repeating the analysis, 
including the t + t~ preselection stage. 

Within the r sample itself the rate of association of different subdetectors to a recon- 
structed charged particle track was studied for different r decay topologies. In candidate 
1-prong decays, the hit association probabilities for the VD, ID and OD were studied for 
tracks which had an associated TPC track segment. The probabilities obtained in data 
and simulation were compared and the relative differences calculated. These differences 
were 0.0% for the VD, +2.3% for the ID and —1.2% for the OD. The equivalent num- 
bers for 3-prong decays were —0.7%, +2.5% and +0.6% respectively. The results for the 
ID were unchanged if in addition the VD was required to be associated to a track, and 
vice-versa. The small fraction of tracks without a TPC track segment showed a level of 
agreement for the proportions of tracks with different subdetectors which was better than 
3% in all cases. 

In simulation, the preselection efficiency was found to be independent of the combi- 
nation of the subdetectors attached to a track at a level below 0.1%. Given the observed 
discrepancies between data and simulation, any effect on the preselection efficiency due 
to this source was O(10~ 5 ). The uncertainties arising from the probability of including 
an ID hit on a track were estimated to be 2.2 x 10~ 5 for both B\ and B% with a full anti- 
correlation and has been included in the tracking systematic uncertainty. The analogous 
uncertainty for the OD was negligible. 

It was observed that the level of Bhabha background was correlated with the existence 
of an OD hit on a track. This was due to bremsstrahlung, and resolution effects in the p ra d 
variable. The level of agreement observed between the data and the simulation implied 



that the Bhabha background was consistent with the estimation made in Section 4.2 



The two-track resolution has been studied by data-simulation comparison of the min- 
imum opening angle in three dimensions between tracks in r decay hemispheres with 
more than one track. Fig. |8|a shows the efficiency in simulation to reconstruct a 3-prong 
r decay as a function of the minimum opening angle. It is flat within 1% except for 
a fall-off of about 5% in reconstruction efficiency for a minimum opening angle below 
3 mrad. As is visible in the distribution of the minimum opening angle in candidate 
3-prong r decays shown in Fig. ||b, less than 1% of 3-prong r decays lie in this region 
and data and simulation are compatible within the statistical precision. A systematic 
uncertainty was attributed to the two-track resolution by varying the efficiency in the 
region with minimum opening angle below 3 mrad by the uncertainty allowed from the 
statistical uncertainty of the data. There was no observable fall-off in efficiency for low 
values of the minimum differences between the azimuthal or polar angles of tracks in a 
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Figure 8: a) The selection efficiency in simulation for candidate 3-prong decays as a 
function of the minimum angle in three dimensions between any two charged particles, 
b ) The distribution of the minimum angle in three dimensions between any two charged 
particles in candidate 3-prong r decays. Dots are data, line is simulation. 



3-prong hemisphere. The minimum opening angle distribution in candidate 5-prong r 
decays is shown in Fig. |9|. 

The association efficiency of the TPC in 3-prong r decays and its dependence on 
the minimum opening angle were studied by comparing the rate of candidate r decays 
containing three tracks with associated TPC hits with the rate of decays containing either 
two tracks with TPC hits plus one track without TPC hits or one track with TPC hits 
plus two tracks without TPC hits. Fig. [H] shows the data over simulation ratio of the 
fraction of candidate 3-prong r decays containing three tracks with associated TPC hits, 
as a function of the minimum opening angle. It shows compatibility with unity in all 
regions except for a small discrepancy of about +3% for minimum opening angles in 
the region of 4 to 7 mrad. The average of the ratio for all minimum opening angles is 
1.0016 ± 0.0008, consistent with unity within two standard deviations. The difference 
from unity can be directly related to the efficiency for a charged particle track to have 
associated TPC hits in 3-prong r decays. The uncertainties on the topological branching 
ratios were estimated by varying the probability in simulated 3-prong r decays to have 
TPC hits associated to a track by 1.6 x 10~ 3 and repeating the analysis. The same 
procedure was carried out for 5-prong decays. 

Fig. |8|b shows some differences between data and simulation in the region between 3 
and 15 mrad, with the simulation tending to lie slightly above the data. This can be due 
to the decay modelling, for which the attendant systematic uncertainties are discussed in 
Section |6.2.4| . However track reconstruction effects cannot be excluded. While systemat- 
ics in the TPC association efficiency are accounted for above, the deviation from unity 
of the ratio shown in Fig [10] can also give an estimate of the rate of track reconstruction 
without the TPC, in particular of VD-ID tracks. An upper estimate of the magnitude of 
this effect for close tracks was derived by integrating up to 15 mrad the deviations from 
unity of this ratio times the number of events in a given bin of minimum opening angle. 
This procedure gave an uncertainty of 1.2 x 10~ 4 on B 3 , together with uncertainties on 
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Figure 9: The minimum angle in three dimensions between any two charged particles in 
candidate 5-prong r decays. Dots are data, line is simulation. 

Bi and B 5 induced via the correlations in the fit. This accounted for potential differences 
between data and simulation in the VD-ID track reconstruction in the small minimum 
opening angle region of 3-prong r decays. 

The reconstructed charge in a r decay hemisphere was studied for different event 
and decay topologies. The charge was not used to classify r decays or select t + t~ 
events, so this provided an indirect cross-check of the track reconstruction. Use was made 
of the constraint that the two r-decay hemispheres arose from r's of opposite charge. 
For events in the data belonging to the 1-1 event class, an estimated (99.88±0.01)% 
of r decays had the charge correctly reconstructed. The rate in the simulation was 
(99.91±0.01)%. For 3-prong decays in data, (98.57±0.09)% had the correctly signed unit 
charge, compared with (98.60±0.03)% in simulation, while (0.54±0.06)% had a charge of 
three with the correct sign as compared with (0.52±0.02)% in simulation. The remainder 
had the wrongly signed charge. In 5-prong decays an estimated (90.5±0.7)% had the 
correctly reconstructed unit charge in data compared with (89.7±0.3)% in simulation. 

6.2.3 Reinteractions and K£ reconstruction 

Uncertainties from the photon conversion reconstruction were estimated by varying 
by their uncertainties the correction factors for the reconstructed and unreconstructed 
conversions, given in Table [l], which were obtained from data test samples as described 
in Section |5.2j . A contribution for the uncertainty in the incident photon rate in the 
dileptonic test samples was included. The resultant uncertainties are dominated by the 
contribution from the unreconstructed conversions. A similar approach was taken for 
the nuclear reinteractions. Here again the uncertainty arising from the unreconstructed 
nuclear reinteractions was the most significant contribution. 

Both the photon conversion and hadronic reinteraction errors contain contributions to 
account for the uncertainty on the rate of fake reconstructions which cause the topology 
to be misidentified. 
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Figure 10: For candidate 3-prong r-decays, the ratio of data to simulation of the fraction 
of candidate decays where three tracks had associated TPC hits, as a function of the 
minimum opening angle. 



The 1-prong r decay class has a number of different subclasses which cover the different 
reconstruction possibilities for an elastic nuclear reinteraction in a 1-prong hadronic r 
decay. In simulation, of the 3.3% of hadrons in r decays undergoing an elastic scatter 
before the TPC sensitive volume 99.4% were correctly attributed to the 1-prong class, 
leading to a classification inefficiency of 2.0 x 10~ 4 for 1-prong hadronic r decays. A 
cross-check of the elastic scatters was performed using a kink algorithm which attempted 
to link up pairs of tracks consistent with both tracks having been produced by a single 
charged particle experiencing a large scatter in detector material. Comparison of the 
rates in data and simulation showed agreement within 20%, which was taken as a relative 
uncertainty on the effect. The effects of fake kink reconstructions by the algorithm had 
a negligible effect. 

The uncertainty attributed to the electron identification had contributions from the 
uncertainties on the efficiencies for conversion rejection and for rejection of 7r° Dalitz 
decays, and from the probability of misidentifying a hadron as an electron. The study 
of the electron identification algorithm and its systematic uncertainties is described in 



Section 5.2 



The (5-ray rate is proportional to the mean atomic number Z of the material, while the 
factors which are applied correct for the number of radiation lengths, which is approxi- 
mately proportional to the mean Z 2 . This gives an ambiguity in the correction for 5-rays 
depending on whether it is due to the wrong mean Z, or to the wrong quantity of mate- 
rial, but with the correct mean Z. The correction factor used in the analysis, Xf /Xf , 
assumes the latter case. The former case would imply a factor of (Xf M / 'Xf M ) l l 2 . This 
was applied and the observed variations on the branching ratios taken as systematic 
uncertainties. 

Only the fraction (~5%) of Kg — ► 7r + 7r~ decays occuring inside the VD, and which 
the Kg reconstruction algorithm failed to identify, gave fake charged primary r decay 
products and hence the incorrect topology assignment. For the region inside the outer 
layer of the VD, the rate of successfully reconstructed Kg 's was 0.91 ± 0.09 times the 
simulation prediction, showing good consistency. The quoted error is purely statistical. 
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By comparison of data and simulation for candidate K§ mesons reconstructed as decaying 
beyond the VD, the reconstruction efficiency and relative rate of photon conversions 
misidentified as K§ were estimated to agree within 15%. The K® reconstruction algorithm 
can reconstruct a fake K§ from two primary r decay particles thus reducing the r decay 
charged multiplicity by two. Simulation studies indicate that the scale of this effect was 
5 x 10~ 5 of the 3-prong rate. The related systematic uncertainty was taken to be equal 
to the size of the effect. 



6.2.4 Exclusive branching ratios and r decay modelling 

The exclusive branching ratios were varied within the uncertainties quoted in the 
Particle Data Listings ||. To take account of any hidden correlations the quoted uncer- 
tainties have been scaled up by a factor of 1.5. The largest single contribution is from 
the decay mode 37r ± 37r i/ r with a high tt° multiplicity and a large relative uncertainty 
on the branching ratio. The uncertainty also included a small contribution due to the 
decay modes r _ — > K~7T~7T + 7T°b> T and r~ — > K~ K + -k^ts: q v t which were not included in 
the simulation. 

The uncertainties associated with the modelling of the 3- and 5-prong decays were 
estimated by correcting the efficiencies taking into account differences between data and 
simulated invariant mass distributions. In addition, the hadronic structure of the 3n 
final state was varied between the default TAUOLA [151 model and that obtained in 



the DELPHI analysis of the 3n structure in r decays p3 |. For the 37T7r° structure the 



parameterisation of Model 1 of [24J] was used and, as a cross-check, the parameterisation 



of 37T7T used in |23| was used to reweight the distributions of the minimum opening angle. 



6.2.5 Trigger, energy scale, r polarisation and simulation statistics 

The trigger efficiency for t + t~ final states was (99.98 ± 0.01)% for events within the 
polar angle acceptance. Studies indicated that the inefficiency was due to events where 
both r's decayed via the r — > {ivv mode JT7[] . It was assumed that the full inefficiency of 
(2 ± 1) x 10~ 4 was contained in the 1-1 event class and the efficiency of this class and of 
the t + t~ preselection were modified accordingly. The associated systematic uncertainty 
was obtained by varying the inefficiency by its error. 

The energy and momentum scales and resolution can affect many of the quantities 
used in the analysis, such as the r + r~ selection variables, invariant masses or the thrust 
(used in the qq rejection). The momentum scale was varied by 0.2%, the electron energy 
by 0.5% and the neutral electromagnetic energy by 0.2%. These variations were obtained 



in a study carried out for the r polarisation measurement |25| . The associated systematic 
uncertainty on the topological branching ratio results was small. 

The 1-prong selection efficiency has a slight sensitivity to the average r polarisation 
because of acceptance effects in r — > i\v T ,Kv T decays due to the p ra d cut used in the 
t + t~ selection. The analysis used the result and uncertainty from the DELPHI analysis 



on t polarisation 25 



The systematic uncertainty due to the limited simulation statistics was also included. 



6.2.6 Other cross-checks and summary 

A cross-check of the fitting procedure was performed by using the simulation itself as 
input to the fit; the results of the fit were compared with the input branching ratios and 
agreement was observed. This was repeated with different input branching ratios, again 
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Source of systematic 


I-prong 


3-prong 


5-prong 


Dilepton background 


110 


-109 


-1 


Cosmic ray background 


5 


-5 


<1 


Four-fermion background 


42 


-41 


-1 


Z — ► qq background 


25 


-24 


-1 


Neural Network qq rejection 


50 


-48 


-5 


Tracking 


157 


-152 


-16 


VD efficiency 


55 


-60 


+6 


Conversions 


126 


-121 


-8 


Inelastic Nucl. reinteractions 


90 


-80 


-10 


Elastic Nucl. reinteractions 


24 


-24 


-2 


Electron identification 


104 


-97 


-7 


5-ray weights 


8 


-8 


<1 


.Kg reconstruction 


5 


-5 


<1 


Exclusive BRs 


228 


-204 


-44 


3-prong decay modelling 


116 


-121 


+ 10 


Trigger 


15 


-15 


<1 


E and p scales 


19 


-20 


+1 


t polarisation 


18 


-19 


+1 


Simulation statistics 


310 


-310 


+31 


Total systematic 


492 


477 


59 


Statistical 


929 


929 


126 



Table 6: Contributions in units of 10' 6 to the systematic uncertainties on B\, B% and B^. 
The uncertainties are signed to show the correlation between the different branching ratios. 
The 1-prong errors are always assumed positive. 

showing good agreement. The reweighting procedure was, where possible, cross-checked 
by direct calculation of the expected variation in the measured branching ratios when 
applying the weights for a given effect. 

In the fit the main contributions to the \ 2 come from the 3-3' and 3'-3' classes. Both 
contain significant background from e + e~ — > qq events. Removing these two classes 
from the fit had an almost negligible effect on the results, within the bounds of expected 
statistical fluctuations. No systematic uncertainty was ascribed to this effect. 

Other cross-checks were performed, including fitting the branching ratios hemisphere 
by hemisphere rather than event by event, using Eqn. For this fit it was practical to 
subdivide the sample even further into hemisphere subclasses dependent on the number of 
VD-ID-TPC(-OD) tracks, VD-ID tracks, number of conversions, nuclear reinteractions, 
etc. The branching ratios obtained with this approach were in excellent agreement with 
those obtained in the event-by-event fit. 

The fits were performed for each year's data separately. Good agreement was found 
for the results in the different years. 

The different systematic uncertainties are summarised in Table |6|. The systematic 
uncertainties are signed so as to give the correlation between the different branching 
ratios. 
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7 Conclusions 

The measurements made of the 1, 3 and 5-prong topological branching ratios were 



B x = (85.316 ± 0.093 sM ± 0.049 
B 3 = (14.569 ± 0.093 stat ± 0.048 
B 5 = (0.115 ±0.013 stat ± 0.006 



syst 



sysj 



70, 



sys, 



70. 



As expected in view of the small value of B$, the branching ratios B\ and -B3 are almost 
fully anti-correlated with a coefficient of —0.98. B 3 and B 5 have a correlation coefficient of 
—0.08, and B\ and B§ also have a correlation coefficient of —0.08. The B\ and B% results 
are consistent with and slightly more precise than the results obtained in the PDG || 
combined fit to all r decay data: B x = (85.32 ± 0.13)%; B 3 = (14.58 ± 0.13)%. The B t 
and B 3 results are more than twice as precise as the existing world averages || of (84.59 ± 
0.33)% and (14.63±0.25)% respectively. The results are in reasonable agreement with, but 
significantly more precise than, the most recent direct measurements, by OPAL [26,27], 
CLEO |28| and ALEPH 



The result on the 5-prong branching ratio is in good agreement with the world aver- 
age of (0.107 ± 0.009)% and the PDG best fit result of (0.099 ± 0.007)% which include 
contributions from OPAL P§, CLEO PJ and ALEPH M measurements. 
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